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Abstract 

It is shown that, for a large class of statistical mixtures, the Wigner-Poisson 
(or Hartree) system can be reduced to an effective Schrodinger-Poisson sys- 
tem, in which the Schrodinger equation contains a new nonlinearity. For 
the case of a zero-temperature one-dimensional electron gas, this additional 
nonlinearity is of the form 4 . In the long wavelength limit, the results ob- 
tained from the effective Schrodinger-Poisson system are in agreement with 
those of the Wigner-Poisson system. The reduced model is further used to 
describe the stationary states of a quantum electron gas and the two-stream 
instability. 



1 Introduction 



Understanding the dynamics of a quantum electron gas is an important issue for a 
variety of physical systems, such as ordinary metals, semiconductors, and even as- 
trophysical systems under extreme conditions (e.g. white dwarfs). Although some 
level of understanding can be achieved by considering independent electrons, a 
more accurate description requires the use of self-consistent models, where electron- 
electron interactions are taken into account. As the treatment of the full iV-body 
problem is clearly out of reach, mean-field models are usually adopted, of which 
the Hartree and Hartree- Fock models are standard examples 0. In the Hartree 
approximation, each electron is described by a one-particle wave function (obeying 
Schrodinger's equation), and the electrostatic force acting on it results from Pois- 
son's equation. (Fock's correction accounts for the parity of the A-particle wave 
function for an ensemble of Fermions, but this correction will not be considered in 
this paper). 

The Hartree model can be written in a more compact, but strictly equivalent, form 
by making use of Wigner functions. The Wigner representation is a useful tool 
to express quantum mechanics in a phase space formalism (for reviews see ||). 
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In this representation, a quantum state (either pure or mixed) is described by a 
Wigner function (i.e. a function of the phase space variables), and the Wigner 
equation provides an evolution equation which is similar to the Vlasov equation, 
well-known from classical plasma physics. We note that, although the Wigner 
distribution satisfies most of the standard properties of probability distributions, 
it cannot be regarded as such, since it may take negative values. The resulting 
self-consistent model is called the Wigner-Poisson (WP) system, and has been 
extensively used in the study of quantum transport [f§-||. 

Despite its considerable interest, the WP formulation presents some intrinsic draw- 
backs : (a) it is a nonlocal, integro-differential system; and (b) its numerical treat- 
ment requires the discretisation of the whole phase space. Moreover, as is often 
the case with kinetic models, the WP system gives more information than one is 
really interested in. For these reasons, it would be useful to obtain an accurate 
reduced model which, though not providing the same detailed information as the 
kinetic WP system, is still able to reproduce the main features of the physical 
system under consideration. 

After discussing the general validity of the WP model, we will derive an effective 
Schrodinger-Poisson (SP) system, which, in an appropriate limit, reproduces the 
results of the kinetic WP formulation. A similar result was recently obtained in 
the mathematical literature [0, |8| , although its physical implications have not been 
fully analysed. In this effective SP model, the Schrodinger equation is nonlinear, 
as it includes an effective potential depending on the modulus of the wavefunction. 
The exact form of this effective potential depends on the specific physical system 
being studied. In order to obtain the effective SP system, we will first derive a 
system of reduced 'fluid' equations by taking moments of the WP system. It will be 
shown that the pressure term appearing in the fluid equations can be decomposed 
into a classical and a quantum part. With some reasonable hyphoteses on the 
pressure term, the fluid system can be closed. Finally, the effective SP system 
will be applied to several physical problems, including linear wave propagation, 
nonlinear stationary solutions, and the two-stream instability. 



2 Coupling parameter for a quantum plasma 

A classical plasma can be said to be collisionless ('ideal') when long-range self- 
consistent interactions (described by the Poisson equation) dominate over short- 
range two-particle interactions (collisions). This happens when the potential en- 
ergy of two electrons separated by an average interparticle distance is small com- 
pared to the average kinetic energy. The potential energy is estimated as E pot = 
e 2 n^ 3 /eo, while the average kinetic energy is simply given by the temperature T 
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(measured in energy units). Here — e is the electron charge, Eq the dielectric con- 
stant in vacuum, and Uq the equilibrium particle density. One defines therefore a 
classical coupling parameter 

-p J^pot e n , 
lc = ~W~ - Z T ' ^ ' 

111 kin e J 

such that the collisionless approximation is valid when < 1. The classical 
coupling parameter can be written in a different way, by introducing the plasma 
frequency, the thermal velocity and the Debye length 

u p = [ ~Er ) ' v t = ( zz ) ' x d = ^ , (2) 
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which are typical inverse time, velocity and length scales for a collisionless plasma. 
With these definitions, the coupling parameter can be expressed as 

which is the inverse of the number of electrons contained in a Debye volume. When 
the condition Tc 1 is satisfied, two-body correlations (collisions) can be ne- 
glected, and the iV-particle Liouville equation can be reduced, via a BBGKY hier- 
archy, to the one-particle Vlasov equation. The Vlasov-Poisson system is therefore 
the standard model to describe classical electrostatic plasmas in the collisionless 
approximation. 

Similarly, it is possible to define a quantum coupling parameter Yq. Let us con- 
sider the case of a completely degenerate electron gas. Now the average kinetic 
energy is given by the Fermi temperature Tp ~ h 2 n^ 3 /m (we neglect irrelevant 
dimensionless constants), so that the quantum coupling parameter becomes 



E pot _ e 2 m ^ 



Tf h 2 e n^ 3 

Notice that, according to Eq. (HI), a quantum electron gas is more ideal at higher 



densities. Using the Fermi velocity Vp = yTp/m, one can define a typical length 
scale for the quantum plasma 

— — — . (5) 

Up 

The quantum coupling parameter can thus be expressed as the inverse of the 
number of electrons contained in a Fermi volume : 

r « 2 = ■ (6) 
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Finally, another expression for the coupling parameter is the following 



1/2 _ fw p 



(7) 




which is valid for any number of dimensions. 

The quantum electron gas is collisionless when Fq <C 1. In this case, the quantum 
A^-body problem can be reduced to a one-particle Wigner equation. The Wigner- 
Poisson system is therefore capable of describing a quantum electrostatic plasma 
in the collisionless approximation. 

The previous results were derived in the limiting cases T ^> Tp (classical) and T <C 
Tp (quantum degenerate). For intermediate temperatures, simple expressions for 
the coupling parameters are not available, but one must expect a smooth transition 
between the two regimes. 

For electrons in metal, we have typically 



n ~ 10 29 m" 3 , v F ~ 10 6 ms" 1 , cu p ~ 10 16 s~\ \ F ~ HT 10 m . (8) 



These values yield a quantum coupling parameter of order unity. Allowing for the 
dimensionless constants we have neglected and the different properties of metals, 
we realise that Tq can be both smaller and larger than unity for typical metallic 



The above values seem to indicate that, as Tq ~ 1, electron-electron (e-e) colli- 
sions cannot be neglected for metals. If that was the case, one should abandon 
one-particle models such as the Wigner or Hartree equations, and resort to the 
full A^-body problem. This is hardly a feasible task. Fortunately, however, the ex- 
clusion principle comes to rescue by reducing the collision rate quite dramatically 
in most cases of interest [|IJ. This occurs when the electron distribution is close 
to the Fermi-Dirac equilibrium at relatively low temperatures. The fundamental 
point is that, when all lower levels are occupied, the exclusion principle forbids 
a vast number of transitions that would otherwise be possible. In particular, at 
strictly zero temperature, all electrons have energies below Tp, and no transition 
is possible, simply because there are no available states for the electrons to occupy. 
At moderate temperatures, only electrons wihin a shell of thickness T about the 
Fermi surface can undergo collisions. The e-e collision rate (inverse of the lifetime 
r ee ) for such electrons is proportional to T/H (this is a form of the uncertainty 
principle, energy x time = const.). The average collision rate is obtained by mul- 
tiplying the previous expression by the fraction of electrons present in the shell of 
thickness T about the Fermi surface, which is ~ T/Tp. One obtains 



electrons. 
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At room temperature, r ee ~ 10~ 10 s, which is much larger than the typical colli- 
sionless time scale r p = u~ l ~ 10~ 16 s. Therefore, for times smaller than r ee the 
effect of e-e collisions can be safely neglected. In addition, it turns out that the 
typical relaxation time scale is r r ~ 1CT 14 s, which is again significantly larger 
than r p . In summary, the ordering 

T p < 7V < T ee , (10) 

implies that a collisionless (Wigner) model is appropriate for relatively short time 
scales. 



3 Derivation of the fluid model 



In one spatial dimension, the Wigner- Poisson system 

df 
dt 



reads as 



df iem 
dx 2nh 



x + 



Xn 



dX dv'e im( - v - v ' )x x 

, / An 



x — 



f(x,v',t) = 



d 2 



dx 2 



~ (/ f dv ~ n o 



(11) 



(12) 



where f(x,v,t) is the Wigner distribution function, 4>(x,t) the electrostatic poten- 
tial, — e and m the electron charge and mass, Eq the vacuum dieletric constant and 
n a background ionic charge. Notice that the one-particle Wigner function used 
here actually represents an iV-particle system. Indeed, the above Wigner-Poisson 
system can be derived from the full iV-body problem via a BBGKY hierarchy, 
neglecting two-body correlations and only keeping the mean Coulomb field |J. 
Further, it is easy to see that, in the limit h — > 0, one recovers the familiar Vlasov- 
Poisson system for classical collisionless plasmas. For simplicity of notation, only 
one-dimensional problems will be treated in the rest of this paper, but the results 
can be readily extended to higher dimensions. 



In order to derive a fluid model, we take moments of Eq. (|TT|) by integrating over 
velocity space. Introducing the standard definitions of density, mean velocity and 
pressure 



n(x, t) = J f dv 
it is obtained 



u(x, t) = — [ fvdv , P(x, t) = m 
n J 



dn d(nu) 

\ — - = , 

dt dx 

du du e d 4> 1 dP 
dt dx md x mn d x 



fv 2 dv 



nu 



(13) 

(14) 
(15) 
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We immediately notice that Eqs. fln|)-(|I5D do not differ from the ordinary evolu- 
tion equations for a classical fluid. This may seem surprising, but in the following 
it will appear that the quantum nature of this system is in fact hidden in the 
pressure term. 

The pressure term may be decomposed into a classical and a quantum part. This 
can be shown as follows. The Wigner distribution for a quantum mixture of states 
ij; a (x,t), each characterized by an occupation probability p a , is written as 



f(x,v,t) 



m 



2nh 



J^Pa J d\ %l)* a ^+^>M i> a (x- 



^imvX/h 



(16) 



where the sum extends over all possible states. The numbers p a , representing 
probabilities, satisfy the relations p a > 0, J2 a Pa — 1- Using the the previous 
expression, one can compute the pressure. After some algebra, one obtains 



P 



4m 



■ 


dlp a 




dx 



%Pa dx 2 



dx 2 



+ 



Aran 



Y,Po 



f , * dip a dip* 



dx 



dx 



1 2 



(17) 



If we represent each state as 

ip a (x,t) = A a (x, t) exp (iS a (x, t)/%) 



(18) 



where A a (amplitude) and S a (phase) are real functions, we obtain P = P° + P Q , 
where the classical P c and quantum P® parts of the pressure are 



2^P a PpA a A 



Iran 




(19) 
(20) 



Notice that P Q only depends on the amplitudes A a , and that for a pure state only 
P® survives. It can be easily shown that P c represents the standard pressure, 
resulting from the dispersion of velocities. To prove this, one has to remember that 
the phases S a are related to the mean velocity u a of each wavefunction through 
the relation rau a = dS a /dx [the u a 's should not be mistaken with the global mean 
velocity u defined in Eq. ([13])]. Thus, by expanding Eq. fllTf ), one obtains after 
some algebra 

^2iaPaA a U a ( ^2 a p a A Q U a \ 



->c 



ran 



n 



n 



(21) 
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With an appropriate definition of averages, we can rewrite the above equation as: 
P c = mn({Ua) — (u a ) 2 ), which is the standard expression for the pressure. The 
contribution P®, on the other hand, is a purely quantum pressure, with no classical 
counterpart. 

In order to close the fluid system, some equation of state, relating P c and P® to 
the density n, must be used. In this paper, we consider a statistical mixture where 
all the amplitudes are equal (but not constant), A a (x) = A(x). This gives, using 
Eqs. (|13|) and (|l6f), n = A 2 . With the help of Eq. ([H]), the quantum pressure 
becomes 



pQ 



2m 



' d_ 
dx 



n 



d 2 



dx 2 



(22) 



For the classical part of the pressure, we make the standard assumption that it 
is some function of the density, P c = P c (n). With these hypotheses, the force 
equation ( |15|) can be written as 



du du 
dt dx 



e d(p 
mdx 



1 OP 



c 



1 dP Q 



Defining the effective potential 

W(n) 

and using the identity 



mn dx mn dx 



dri dP c {n') 



n 



dn' 



1 8P Q h 2 3 (d 2 {^n)jdx 2 



mn dx 



2m 2 dx 



(23) 



(24) 



(25) 



the force equation ( 23 ) reduces to 

edd> 1 dW 



du du 
dt dx 



+ 



h 2 d (d 2 {jn)jdx 2 



mdx mdx 2m 2 dx 



n 



(26) 



Now comes the crucial point : it is possible to combine Eqs. ([14]) and (|26|) into an 
effective Schrodinger equation. Indeed, let us define the effective wavefunction 

* = yjn(x,t)exp (iS(x, t) /%) , (27) 

with S(x,t) defined according to mu(x,t) = d S(x,t)/d x. We obtain that *&(x, t) 
satisfies the equation 

<9\& h 2 d 2 ^/ 

!fi a? = -ra-^* + ^' (28 > 



7 



This is a nonlinear Schrodinger equation, as the effective potential W depends on 
the wavefunction through Eq. (pi|) , where n = \^\ 2 . Separating Eq. (p8|) into 
its real and imaginary parts, we indeed find the continuity (13) and force (p6|) 
equations. Finally, the complete effective SP system is composed of Eq. ( p8l) and 
the Poisson equation 

d 2 6 



-(|vl/| 2 -n ) . (29) 



dx 2 e 

To summarize what we have achieved so far, we notice that, in general, the dynam- 
ics of a statistical mixture must be treated with the full Wigner-Poisson system, 
or, equivalently, with a set of Schrodinger equations, coupled by Poisson's equa- 
tion (Hartree's model). In the present section, we have shown that one can reduce 
the problem of quantum transport to a single nonlinear Schrodinger equation plus 
Poisson's equation. Also, notice that the nature of the interaction (electrostatic in 
our case) is not of essential importance. The main result is that we can reduce the 
(phase space) Wigner equation to a (real space) nonlinear Schrodinger equation. 

The two hypotheses used for this reduction are : (a) all states composing the mix- 
ture have the same amplitude [which leads to Eq. ( P2| ) for the quantum pressure]; 
and (b) the equation of state for the classical pressure is P c = P c (n). Hypothesis 
(b) is the standard fluid closure, and needs no further comment. Hypothesis (a) 
means that all electrons are distributed in space according to the same probability 
distribution n(x) = A 2 (x). What distinguishes the electrons from one another is 
their phase S a , and therefore their velocity u a . This approximation can be viewed 
as a first step beyond the standard homogeneous equilibrium of a Fermion gas, for 
which each state is represented by a plane wave 

if) a (x,t) = Aexp (imu a x/h) , 

with the amplitude A and the velocity u a spatially constant. In our approximation, 
both the amplitude and the velocity can be spatially modulated, although we still 
restrict ourselves to the case where the amplitude is the same for all states. This 
appears to be a reasonable closure assumption for systems that are not too far 
from the Fermi-Dirac equilibrium. 



4 Applications 

As a relevant example of the above theory, we consider a zero-temperature one- 
dimensional electron gas, with Fermi velocity Vp and equilibrium density uq. In 
this case, the classical pressure is 
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(Notice that the term 'classical' is somewhat inappropriate here, as P° will contain 
Planck's constant through the Fermi velocity). We also note that the Fermi velocity 
in one spatial dimension 

V F = (31) 



m 



is proportional to n 0l whereas in three dimensions vp oc n^ 3 . 

Using Eq. (|30|) , the effective potential defined in Eq. (|4]) turns out to be 



W=?%W. (32) 

Notice that the effective potential is repulsive, and tends to flatten the electron 
density. This is quite natural, as W derives from the pressure P c , which in turn 
is a manifestation of the dispersion of velocities in a Fermion gas. When the gas 
is at equilibrium, W ~ n ~ const., and this term has no effect. 

We also point out that a similar nonlinear Schrodinger equation with a |^/| 4 - 
dependent potential has recently been derived in the study of low-dimensional 
Bose condensates [fLOR . We stress, however, that such a Boson- Fermion duality 
only applies to one-dimensional systems. For a D-dimensional Fermion system, 
the classical part of the pressure has the form P ~ n( D +W D , so that the effective 



potential becomes W ~ n 2 ^ D BTT 



4.1 Linear wave propagation 



As a first application, let us study wave propagation for the effective SP system 
(p8|)-(p9|) with W given by Eq. fl32|). Linearizing around the homogeneous equi- 
librium ^ = y/nQ, e<p = mvp/2, we obtain the following dispersion relation (for 
waves with frequency uj and wave number k) 

u* = ul + k\\ + ^ . (33) 

For vf = 0, we recover the dispersion relation of the standard SP system |J. 
Equation (|33| ) can be written in dimensionless units by using Eqs. (||) and (0), 
which are valid both in one and three spatial dimensions 

- 2 =l + e\l+ k -^Y Q . (34) 

Note that quantum mechanical effects (dispersion of the wave packet) are first 
order in the coupling parameter Tq, whereas quantum statistical effects (Fermi- 
Dirac distribution) appear at leading order. 
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We want to compare this dispersion relation to the one obtained from the complete 
WP system (p])-(12), which, in the most general case, reads as [||, [L^ 



ujI f f Q (v) dv , , 

n J (uj - kvf - ft 2 fc 4 /4m 2 ' 1 ; 

In our case, fo(v) is given by the Fermi-Dirac distribution for a zero-temperature 
one-dimensional electron gas at equilibrium, i.e. fo(v) = n /2v F if \v\ < v F and 
fo(v) = if \v | > Vp. Substituting into Eq. fl35"D , one obtains (without any further 
approximation) 

— 2 =^ coth — + k 2 \ 2 F + — ^ T Q , (36) 



CO 2 UJ 2 \ CO 2 

p p \ p 



where 



n 2 _ hk 3 v F 

u 2 muj 2 



k'X% . (37) 




Now we expand the first term on the right-hand side of Eq. (^) in the long 
wavelength (fluid) limit Q u p . Using the expansion a;coth(x) = 1 + x 2 /3 — 
x 4 /45 + . . ., one obtains 

^ + ^)r -^ 12 A!?rJ + .... (38) 

This is a double expansion in powers of the parameters Tq and k\ F . The colli- 
sionless regime is in principle characterized by Tq 1, although, as was seen in 
Sec. 2, electron-electron interactions can be neglected even when ~ 1, as is 
the case for metals. On the other hand, the fluid regime is characterized by small 
wave numbers (Q <^ uj p ). Indeed, keeping terms to fourth order in k\ F , Eq. (j3S|) 
reduces to the dispersion relation for the effective SP system, Eq. (|3"4]). This is 
a further indication that the effective SP system is a good approximation to the 
complete WP system for long wavelengths. 

We also note that for — >• 0, the dispersion relation reduces to 

uo 2 = uo 2 p + k 2 v 2 F . (39) 

This is exactly the dispersion relation obtained from the classical Vlasov-Poisson 
system with a zero-temperature Fermi-Dirac equilibrium. In other words, when the 
quantum coupling parameter is vanishingly small, a classical dynamical equation 
can be used, as the only quantum effects come from the Fermi-Dirac statistics. 
This situation may apply to extremely dense astrophysical systems such as white 
dwarfs. 



10 



4.2 Stationary solutions 



As a second illustration, we use the present formalism to describe the stationary 
states of the electron gas [^3[ . This result is more easily obtained by using the 
fluid version of our model. In the time-independent case, the continuity equation 
(|HD and the force equation possess the following first integrals 

,o ^ mu 2 , ft 2 <i 2 y4 . , 

J = AV B =_-* + H r___, (40) 

where A = y^n. The first integrals in Eq. (|40|) corresponds to current (J) and 
energy (E) conservation. We can always choose E = by a shift in the electrostatic 
potential. In this way, we can reduce the description of the stationary states to a 
set of second-order nonlinear ordinary differential equations for the amplitude A 
and the electrostatic potential <fi. For a zero-temperature one- dimensional electron 
gas, the effective potential W is given by Eq. (p2|); thus from Eqs. (pEDP , fl29|) we 
get 

~d 2 A (mJ 2 . , mv% tC \ , iS 

R V = m {iv- 2eA * + if A )' < 41 » 

% - 7y- n " ] - (42) 

Notice that, if the amplitude A(x) is a slowly varying function of x, the second 
derivative on the left hand-side of Eq. can be neglected. With this assumption, 
Eq. ( f4l"D reduces to an algebraic equation, which can be solved for A, and the 
result plugged into Eq. (^2]). This becomes a nonlinear differential equation for 
the electrostatic potential, which is nothing but the Thomas- Fermi approximation 
to our model. 

It can be easily verified that the J = case cannot sustain small-amplitude, 
periodic solutions. Hence, we assume J = hqUq with uq ^ and introduce the 
following rescaling 

uj„x j A ecj) 
x = — — , A = 



u y/n^ ' mul 

H = V F = ^. (43) 

We obtain, in the transformed variables (omitting the circumflex for simplicity of 
notation), 



= + (44) 



d 2 ct) 
dx 2 



A 2 -I, (45) 
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a system that only depends on the rescaled parameters H and Vp. Note that the 
quantum coupling parameter can be written as Tq = H/Vp. 

It is interesting to perform a linear stability analysis in order to see in what con- 
ditions the system supports small amplitude spatially periodic solutions. Writing 

A = 1 + A'exp(ikx) , <p = (1 + V|)/2 + <p' exp(ikx) , (46) 

and retaining in Eqs. (f|4|)-(f|g) only terms up to first order in the primed variables, 
we obtain the relation 

H 2 k 4 -4(1 - V^)k 2 + 4 = 0. (47) 
This second degree equation has solutions 



„ 2(1 - V 2 ) ± 2 J (I - VpT - H 2 
k 2 = - F - ^ - . (48) 

Clearly, spatially oscillating solutions only exist when k 2 is real and positive, which 
yields the condition 

V 2 < 1 - H , (49) 

or equivalently 

m«Q > mvp + hujp . (50) 

This expression sets a lower bound on the speed uo, below which no oscillating 
stationary solution can exist. 



4.3 Two-stream instability 

A classical plasma composed of two counterstreaming electronic populations with 
velocities ±u can give rise, for certain wave numbers, to an instability. In a 



previous paper | 13| , we have shown that quantum effects modify the dispersion 
relation, and give rise to a new instability branch. These results were obtained 
by neglecting the effects of quantum statistics, and are therefore valid in the limit 
Vp Uq- Here, we perform the same calculations for finite values of Vf. 

We consider two electronic populations, which are both distributed according to a 
zero-temperature Fermi-Dirac equilibrium, but with average velocities ±n . The 
motionless ions provide a neutralizing background. The dispersion relation for 
such a two-stream plasma can be found in the following way. For a single stream 
propagating at velocity ±w , our fluid model yields the following dielectric constant 
(thus valid for long wavelengths) 

uj 2 

e ± (k, u) = 1 - ( uTkuQ y_ k 2 P v 2_tf k i/ 4m 2 ■ ( 51 ) 
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Setting e±(k,uj) = leads to the dispersion relation found previously, Eq. (|33|), 
with the appropriate Doppler shift. The dielectric constant for the two-stream case 
is found by averaging the contributions from each stream e(k, uo) = (e+ + e_)/2. 
Using the normalization of Eqs. (fE|), we obtain 



e{k, uj) = 1 



1/2 



1/2 



(cu + k) 2 - k 2 Vp - H 2 k 4 /4 (u - k) 2 - k 2 Vj - H 2 k 4 /A ' 



(52) 



Setting e{k,uS) = 0, we obtain the dispersion relation for the two-stream plasma 

u 4 - (l + 2k 2 (l + V 2 ) + ^)u; 2 



t 2 /-i 1,2 H 2 k 2 \ ( t r2\ 1 2 H 2 k 4 

k 1 - Vi - 1,1 ' 1 ■ - w - 



4 



1- (1- V^)fc 2 + 



4 



0. 



(53) 



Notice that for Vf = we recover the dispersion relation obtained in ||13|| . Solving 
for uj 2 , one obtains 



= - + k 2 1 + vl H ± - 

2 V ^ 4/2 



1 + 8fc 2 ( 1 + 2k 2 Vl + 



if 2 A; 



1/2 



(54) 



The solution for c<j 2 has two branches, one of which is always positive and gives 
stable oscillations. The other solution is negative (u 2 < 0) provided that 



[H 2 k 2 - 4(1 - V$)][H 2 k 4 - 4(1 - Vl)k 2 + 4] < . (55) 

We immediately notice that, if Vp > 1, Eq. (|55|) is never verified, and therefore 
there is no instability. This is a quite natural result. Indeed, mathematically, the 
instability is due to the fact that the two-stream velocity distribution has a 'hole' 
around v = 0. When Vp > 1, the hole is filled up, and no instability can occur. To 
put it differently, there can be instability only when the equilibrium distribution 
is a non-monotonic function of the energy, which ceases to be true when Vp > 1. 

When Vp < 1, Eq. ( |55D bifurcates for H — 1 — Vp. If H > 1 — Vp, the second factor 
is always positive, and instability occurs for H 2 k 2 < 4(1 — Vp). If H < 1 — Vp, 
there is instability if either 



or 



< H 2 k 2 < 2(1 - V£) - 2^(1 - V$) 2 - H 2 , 

2(1 - V 2 ) + 2^(1 - Vl) 2 - H 2 < H 2 k 2 < 4(1 - V|) 

13 



(56) 
(57) 



This yields the stability diagram plotted on Fig. 1, which is similar to the one 
obtained in the limiting case Vp = 0. The presence of a finite Fermi velocity has 
the effect of reducing the region of instability. Numerical simulations yield similar 
reult to those observed in the Vf = case, which are reported in [|l3f| . 



5 Conclusion 



In this paper, we have first established the conditions of validity of the Wigner- 
Poisson system. Subsequently, by taking moments of the Wigner equation, we have 
derived an effective Schrodinger- Poisson system that captures the essential features 
of a quantum electron gas. In the long wavelength limit, this model correctly re- 
produces the results of the linear analysis of the Wigner-Poisson system. The 
advantages of the effective SP model are manifold: it is local in space (compared 
to the nonlocal WP system); it is cast into the ordinary space, rather than the 
phase space; and it has a straightforward interpretation in terms of fluid quantities. 
Furthermore, it is easily amenable to numerical studies, given the abundance of 
accurate numerical techniques for the Schrodinger equation (in comparison, numer- 
ical methods for the Wigner equation |0| are much scarcer and more cumbersome 



to implement). The crucial points in the derivation of the model are : (a) the 
decomposition of the pressure into a classical and a quantum contribution, and (b) 
the restriction to an appropriate class of statistical mixtures (composed of states 
with same amplitude, but different phases). We believe that this class is wide 
enough to describe a significant range of relevant physical systems. 

For the case of a completely degenerate electron gas, the effective SP model can 
be put in a particularly simple form, in which the Schrodinger equation exhibits 
a l^l 4 nonlinearity. This model has been applied to the study of linear wave 
propagation, nonlinear stationary solutions and the two-stream instability. The 
simplicity of the resulting system of equations makes it a useful tool for the study 
of quantum transport in solid state plasmas. 
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Figure Captions 

Figure 1 : Stability diagram for the two-stream plasma, with Vp = 0.7 (solid 
lines) and Vp — (dashed lines). The curves correspond to Eqs. (^)-(|57|). For 
both cases, the region of the plane containing the H 2 axis is unstable. 
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